Single-Chain Magnets from sharp to broad domain walls 
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We discuss time-quantified Monte-Carlo simulations on classical spin chains with uniaxial 
anisotropy in relation to static calculations. Depending on the thickness of domain walls, controlled 
by the relative strength of the exchange and magnetic anisotropy energy, we found two distinct 
regimes in which both the static and dynamic behavior are different. For broad domain walls, the 
interplay between localized excitations and spin waves turns out to be crucial at finite temperature. 
As a consequence, a different protocol should be followed in the experimental characterization of 
slow-relaxing spin chains with broad domain walls with respect to the usual Ising limit. 
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I. INTRODUCTION 

The interest in the physics of domain walls (DWs) 
in Id magnetic systems has been renewed by the capa- 
bility of controlling their motion by means of an elec- 
tric currenti. The technological relevance of this topic 
mainly derives from the possibility of employing DWs in 
novel magneto-storage and spintronic devices^. From a 
more fundamental point of view, the synthesis of the first 
slow-relaxing spin chains^r™ gave the chance to reconsider 
thermally-induced DW diffusion in Id classical spin mod- 
els. The systems displaying such a behavior have been 
named Single-Chain Magnets (SCMs), by analogy with 
Single- Molecule Magnets (SMMs)^: both classes of mate- 
rials show a magnetic hysteresis at finite temperature due 
to slow dynamics rather than cooperative 3d ordering. 
Thanks to the remarkable property of being bistable at a 
molecular level, SMMs and SCMs have been proposed as 
possible magnetic storage units. However, the advantage 
of having identical units whose arrangement might - in 
principle - be tailored with chemical methods is coun- 
terbalanced by the relatively poor thermal stability. In 
fact, the relaxation time becomes macroscopic only at 
temperatures of the order of one Kelvin or lower. The 
quest to improve thermal stability of molecular magnets 
has called for a better understanding of their physical 
properties. The investigation and development of SCMs 
has become an independent and very active field of re- 
search during the last decadeS"— . At odds with SMMs, 
which can be considered Od magnetic systems, SCMs de- 
velop short-range order over a distance comparable to 
the correlation length, ^. The latter is expected to di- 
verge when the temperature approaches zero. However, 
defects and lattice dislocations - whose role is particu- 
larly dramatic in Id - typically hinder the divergence of 
^ below a certain temperature. The physics of SCMs 
is thus different depending on whether the correlation 
length exceeds or not the average size of connected spin 
centers. Finite-size effects will be neglected in our theo- 
retical investigation. This means that our results should 



apply to the temperature region for which the correla- 
tion length is smaller than the average distance between 
two defects in a real spin chain. Under these hypothe- 
ses, some DWs will be always present in the system at 
finite temperature. A simple random- walk argument re- 
lates the relaxation time, r, to the correlation lengtbii: 
within a time r a domain wall performs a random walk 
over a distance ^. The characterization of SCMs basically 
consists in measuring both these quantities and r) as 
function of temperature. In ideal cases, the observed be- 
havior is then reproduced by fitting the parameters of an 
appropriate spin model to the experimental data. 

Following the experimental procedure, we studied 
the temperature dependence of the correlation length 
and the relaxation time in a representative model for 
classical spin chains with uniaxial anisotropy. ^ has 
been computed with the transfer-matrix technique^^^^-i . 
The relaxation of the magnetization and DW diffu- 
sion have been studied by using time-quantified Monte 
Carlo (TQMC)i^"ii. The latter is a recently devel- 
oped algorithm which, in classical spin systems, simu- 
lates a dynamics equivalent to the stochastic Landau- 
Lifshitz-Gilbert equation. TQMC has been previously 
employed to model classical spin chains in the high damp- 
ing limiti^ii^, but the most recent developments allow 
using it for low-damping calculationsiS. 

Most of SCMs show a large single-site magnetic 
anisotropy. The reference model is thus represented by 
the Id Ising model or - more specifically - its kinetic ver- 
sion proposed by Glauber—. In order to better capture 
the physics of experimental systems, the Glauber model 
have been extended to take into account finite-size ef- 
fects^"—, ferrimagnetism^l^ the effect of a strong exter- 
nal magnetic field^^ and the reciprocal non-collinearity of 
local anisotropy axe a'^'^^i'^" . Here we report about a lin- 
ear ferromagnetic spin chain and we address the question 
of how physics changes when the single-site anisotropy is 
progressively reduced. We found that the Glauber sce- 
nario needs to be revisited for SCMs which do not possess 
a large single-site anisotropy so that DWs acquire a finite 
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thickness (more than one lattice unit). Our theoretical 
predictions are in agreement with the few available exper- 
imental results on molecular spin chains with low single- 
site anisotropy^^'"^^. Note that in metallic nonowires (Co, 
Ni, Fe, Permalloy) DWs always extend over several lat- 
tice units. Therefore, for materials traditionally used 
in magneto-storage manufactory the Glauber's picture 
of magnetization reversal and relaxation is not expected 
to hold exactly, nor the correlation length is expected 
to show an Ising-like behavior over a wide temperature 
range. 

In section [TTl we present the model and define the 
physical questions that we want to address. In section 
mil we study the temperature dependence of the corre- 
lation length by means of the transfer-matrix technique 
and Polyakov renormalization. In section IIVI we intro- 
duce the time-quantified Monte Carlo method and study 
the temperature dependence of the relaxation time in the 
broad- wall regime. In section |Vl we use TQMC to study 
temperature-induced diffusion of both sharp and broad 
DWs. In section IVll we provide some phenomenological 
arguments which support our findings and discuss how 
they compare with experiments. In the conclusions, we 
summarize our main results and highlight their relevance 
for further work which could possibly include an electric 
current driving DW motion. 
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FIG. 1. Domain-wall energy (in J units) as a function of 
the ratio D/J, adapted from Ref. [3^ : discrete-lattice cal- 
culation Edw (solid line), continuum- limit solution 2\/2DJ 
(dashed line). Crosses and circles correspond to the ratios 
used in transfer-matrix and time-quantified Monte-Carlo cal- 
culations respectively. The vertical line indicates the value 
{D/J = 2/3) at which the transition from broad- to sharp- 
wall regime occurs. Insets: sketch of the spatial dependence 
of Si as a function of i in the broad-wall (up left) and in the 
sharp-wall (down right) regime. 



II. THE SYSTEM 

As a reference model for SCMs, we consider the fol- 
lowing classical Heisenberg Hamiltonian: 
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where D represents the anisotropy energy, J the exchange 
coupling and H an external applied field. Each spin vari- 
able Si is a three-component unit vector associated with 
the i-th node of the lattice. In this paper either periodic 
or open boundary conditions will be considered depend- 
ing on the calculation. D and J are assumed positive so 
that Hamiltonian ([T]) describes a spin chain with uniaxial 
anisotropy pointing in the z direction. This model and 
variations of it have been investigated extensively from 
the theoretical point of viewi^i^"— . Moreover, Hamilto- 
nian ([T]) has been employed to reproduce the experimen- 
tal behavior of some SCMsi^i^^. Many physical proper- 
ties are related to the energy increase due to the creation 
of a domain wall (DW) in one of the two ground-state, 
uniform configurations with S^' =±1 (Vi). In Fig. [1] we 
plot this energy obtained from a discrete-lattice calcula- 
tion Edw (solid line) for different values of the ratio D/J. 
For D/J < 2/3, the minimal DW energy is realized by a 
spin profile in which several spins are not aligned along 
the easy axis, z (see the sketch in the inset up on the 
left). In this case, DWs spread over more than one lattice 
spacing: broad DWs. For £>/ J > 2/3, the minimum DW 
energy is realized if the transition between — +1 to 
5,f = —1 occurs within one lattice spacing; in this way all 
the spins are aligned along the easy axis: sharp DWs (see 
the sketch in the inset down on the right). The transi- 
tion between broad- and sharp- wall regime is highlighted 
by a singularity in the log-linear plot, which evidences 
a different functional dependence of Sdw on the param- 
eters D and J for anisotropy-to-exchange ratios smaller 
or larger than 2/3. For D/J-^1, the analytic expression 
2V2DJ can be obtained by minimizing the DW energy 
in the continuum- limit approximation^^ (see Eq. 
This function is plotted in Fig. [1] as a dashed line. Ap- 
proaching the transition ratio, D/J = 2/3, from below 
the continuum- limit prediction starts to deviate from the 
discrete-lattice calculation, £dw (solid line). This is rea- 
sonable since the continuum approximation is more ac- 
curate for smaller anisotropy-to-exchange ratios, namely 
the broader DWs are. For D/J> 2/3, the DW energy 
takes the constant value Edw = 2 J. The transition region 
where none of the two analytic expressions holds exactly 
is narrow, meaning that the function 
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D/J < 2/3 
D/J > 2/3 



(2) 



describes the DW energy accurately for most values of 
D and J. In Fig. [U symbols represent the DW energy 
Edw for the values of D/J at which static (circles) and 
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dynamic (crosses) calculations have been performed. 
The main question we want to address is the following: 
are SCMs ruled by different laws depending on whether 
DWs are sharp or broad? Our results show that physics 
is significantly different in these two limits. This fact 
will eventually affect the experimental characterization 
of SCMs which can be modeled by Hamiltonian ^ (or 
variations of itii^iiSi^-32 ) . 

Two key quantities characterizing a specific SCM are the 
correlation length, ^, and the relaxation time, r. For our 
study, the most relevant correlations are those relating 
to spin projections along the easy axis. Therefore, the 
correlation length shall be defined from pair-spin corre- 
lations along z as: 



(3) 



(where (...) stands for thermal average). The static sus- 
ceptibility along z scales with the correlation length as 
follows: 
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The relaxation time can be obtained from the dynamic 
susceptibility. 



X(^, T) - 



X{T) 
1 — iujT ' 



(5) 



where uj is the frequency of the oscillating applied field 
and x(r) is the static susceptibility. Both the real and 
the imaginary part of x('^i display a maximum for 
LJT — 1. In the following we will use an alternative defini- 
tion of T in terms of the relaxation of the magnetization. 
Eqs. (jl]) and ([5]) relate ^ and r to measurable quanti- 
ties, the static and dynamic susceptibility. Based on a 
random- walk argument the relation between the cor- 
relation length and the relaxation time 
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is usually assumed^"^*' for such temperatures that £^<N 
(bulk regime). Dg is the DW diffusion coefficient but it 
can also be interpreted as the attempt frequency for a 
single-spin flip. The temperature dependence of Ds will 
be discussed in more details in Sect. |Vl Within the ki- 
netic Ising model, it is possible to deduce Eq. Q analyt- 
icallyiii2£. The basic experimental characterization of a 
SCM essentially reduces to determining the temperature 
dependence of the three quantities involved in Eq. t, 
^ and Dg. The scenario is well-established in the sharp- 
wall regime, in which one expects that all these quantities 
obey a thermally activated mechanism: 
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fcs = 1 will be assumed throughout the manuscript. 
Eq. ^ relates ^ and t - two experimentally accessible 
quantities - with each other so that the following relation 
between energy scales holds: 
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(8) 



The above relation has been confirmed by several exper- 
imental works on sharp- wall SCMs^"— using the reason- 
able assumption Aa ~ D proposed in Ref. and H^. 
In Sect. |V] we will see that our numerical results confirm 
the validity of Eq. ([H]) for sharp DWs. On the contrary, 
the scenario turns out to be significantly different in the 
broad-wall regime. In fact, for £)/J<2/3, we found that: 

• A^ is not constant but is reduced by more than 30% 
of the value that it takes at low T with increasing 
temperature. 

• Ds does not follow a thermally activated mecha- 
nism as suggested by Eq. ([7]). 

These two important findings indeed affect the experi- 
mental characterization of a SCM for which D/J<2/3. 



III. STATIC PROPERTIES 

A. Transfer matrix calculations 

The static properties of a classical spin chain can 
be computed efficiently with the transfer-matrix (TM) 
methodi3A40. Here we use this technique to compute the 
thermodynamic properties of an infinite chain {N — > oo) 
but finite systems could also be considered'^'*. The essen- 
tial ideas of the TM approach are recalled in Appendix El 
In particular, we used this method to compute the cor- 
relation length. The dependence on r can be eliminated 
from Eq. ^ by summing over all the lattice separations 
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In practice, we evaluated the summation numerically and 
inverted the previous formula as follows 
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Fig. [5] shows the logarithm of the correlation length vs. 
the energy of the DW divided by the temperature as ob- 
tained by TM calculations. Different curves correspond 
to different values oi D/J. Since we set J= 1, the cases 
with D = 1 — 4 fall in the sharp-wall regime, while the 
cases with D — 0.1 — 0.6 correspond to the broad- wall 
regime (see Fig. [T|). We will discuss the behavior of ^ in 
these two regimes separately. 

Sharp-wall regime - In this regime, the DW energy 
amounts to — 2 J. This value has been used to plot 
the logarithm of the correlation length as a function of 
Ew/T for £) = 1 — 4 in Fig. ^ At low temperatures, all 
the solid curves have the same slope, equal to one, re- 
vealing that In(^) ~ Ew/T. A reference line with slope 
one is plotted with short dashes. For D = l, some devia- 
tions from this straight line occur at high temperatures. 
Later on, we will show that in the broad-wall regime the 
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FIG. 2. Color online. Logarithm of the correlation length 
obtained by TM calculations as function of the inverse tem- 
perature in units of the DW energy Ew (defined in Eq. 0) for 
J = 1 and different values of D: Dotted lines, broad domain 
walls, L» = 0.1 (black), 0.2 (red), 0.3 (green), 0.4 (blue), 0.5 
(cyan), 0.6 (magenta). Solid lines, sharp domain walls, D=l 
(black), 2 (red), 3 (green), 4 (blue). Dashed horizontal lines 
give the reference for typical lengths of realistic spin chains. 



interplay between spin waves and DWs leads to an effec- 
tive decrease of with increasing T. It is reasonable 
to think that a similar phenomenon may take place in 
the sharp-wall regime as well, when the transition ratio 
D/J = 2/3 is approached from above. The horizontal 
dotted lines indicate the values ^ = 10^, 10'^, 10^. Due 
to lattice dislocations or impurities (in molecular com- 
pounds) or intrinsic problems in the deposition procedure 
(in mono-atomic nanowires^), the average length of spin 
chains is typically of 10^ — 10'' magnetic centers in real 
SCMs. The actual length of a spin chain sets an upper 
bound to the low-temperature divergence of ^. Such an 
upper bound can be reduced by introducing additional 
non- magnetic impurities^^. Most of the SCMs reported 
in the literature fall in the sharp-wall regime^"—. For 
them, in the temperature range where the correlation- 
length divergence is relevant (^3>1), can be assumed 
equal to 2 J, independently of T. 

Broad-wall regime - For I^/J < 2/3, the DW energy is 
described well by the function £yj~2^/2DJ (see Eq. ([2|)). 
As already pointed out when discussing Fig. [TJ this ex- 
pression obtained in the continuum limit deviates from 
the discrete-lattice calculation in the vicinity of the tran- 
sition region from broad to sharp DWs. In order to 
make the comparison with experiments easier, in the 
main frame of Fig. [2] the correlation length is plotted 
vs. Ew/T, with £yj = 2\/2DJ for all values of D falling 
in the broad- wall regime (instead of using Sd-un value 
obtained from the discrete-lattice calculation). At low 
temperature the slope of all the dotted li nes is about 0.9, 
which suggests an effective A^ = 1.8\/2_D J. This ten 
per cent of reduction with respect to the DW energy can 
be accounted for using the low-temperature expansion of 



the correlation given in Ref. '36': ^ ^ {T /sw) exp {ew/T). 
More strikingly, all the dotted curves bend when the 
temperature increases till their slope becomes roughly 
0.6 at high temperature. A tentative fit of the temper- 
ature dependence of ^ using the corresponding formula 
in Eq. ([7]) would give a value of the activation energy 
1.2^/2DJ < A^ < 1.8V2DJ depending on the tempera- 
ture range in which the fitting has been performed. This 
is nothing but the standard procedure followed in the ex- 
perimental characterization of SCMe^"^°. Normally, the 
extracted energy is then related to the barrier observed 
in the relaxation time by means of Eq. ([5]). It is worth 
remarking that the temperature ranges where A^ and At- 
are extracted usually do not overlap. In fact, to measure 
f the characteristic time scale of the experiment has to be 
much longer than the relaxation time; while for measur- 
ing T the characteristic time scale of the experiment has 
to be comparable to the relaxation time itself, or shorter. 
Thus, the effective variation of A^ with T - which may 
reach 30% of its value - has to be considered when one 
tries to relate this quantity to A^ by means of Eq. (|5]). In 
particular, we remark again that A^ can only be accessed 
at relatively high temperature because finite-size effects 
or three-dimensional inter-chain interactions prevent the 
correlation length from diverging indefinitely'^^'-^. 
As an example of how this first theoretical result may 
be used to better characterize real SCMs, we refer to the 
system reported in Ref. Isil . That spin chain is better 
described by the Seiden modeP^ with anisotropy rather 
than Hamiltonian ([1]). Apart from a multiplicative factor 
two in front of the exchange-energy term, the Hamilto- 
nian of the Seiden model and the one we study here take 
the same form in the continuum limit ''^ (see Eq. (|lip ). 
Yet, the experimental system shows A^^^ ~ 60 K and 
J = 108.5 K (adapting the fitted value to our model). 
Thus 

1. assuming that the fit of A|^^ has been performed 
in the region where A^ — O.Gsw, one would obtain 
e^~100 K and, accordingly, _D=:11.5 K; 



2. assuming that in the fitted region A^ =0.9et, 
would get D = 5.1 K. 



one 



The first estimate of D falls in the suggested range 
D — 8.2 — 17.2 K (estimated by EPR measurements 
on equivalent isolated magnetic units^) while the sec- 
ond estimate is clearly wrong. The experimental cor- 
relation length displays a clear exponential divergence 
in the range 25 K < T < 50 K, which corresponds to 
2 < e^i/T < 4 {ii — 100 K is assumed). Indeed, these 
values of the reduced variable e^/T relate to the high- 
temperature region where A^~0.6£u, for D/J ~ 0.1 (see 
Fig.©. 

In the following we will give a justification of the im- 
portant variation of A^ with T observed in the broad- wall 
regime in terms of Polyakov renormalizationiS. 
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B. Polyakov renormalization 

We focus now in the broad-wall regime in which 
decreases with increasing temperature. A deeper insight 
in such a phenomenon can be achieved by considering 
the effect of spin-wave renormalization on the coupling 
constants D and J. To this aim, we work with the con- 
tinuum version of the Hamiltonian in Eq. ([l}: 



n 



-JN 



i\dJ\^-DiS%x)f 



dx , 



(11) 



where a unitary lattice spacing has been assumed. Fol- 
lowing Polyako v^^i*^^ , we represent S{x) as a superposi- 
tion of a field fluctuating over short spatial scales, 4>{x), 
and a field varying smoothly and over large spatial scales, 
n(x). More explicitly, we write 



S{x) = ri{x) Vl-(/)2(x) + (t>{x). 



(12) 



Requiring 15(2;) | = 1 and |Ti(a;)| = 1, one necessarily has 
n{x) ■ <j){x) = 0; thus can be expressed on a basis 

orthonormal to n{x): 



with a = 1, 2 



(13) 



with |ea(a;)| = 1. The terms appearing in Hamilto- 
nian are affected by the averaging procedure over 
the Six) field as follows: 



=[l-i{<Pl)]{n^f + {^l). 



(14) 



The latter equations represent a well-known result^'^, 
however we will recall their derivation in Appendix [B] for 
convenience. 

In the isotropic Heisenberg chain {D~{)) and in the easy- 
plane case (_D < 0) excited spin waves suffice to destroy 
the long-range order present in the ground state at finite 
temperatures (because of the existence of a Goldstone 
mode). When D > 0, the spectrum of spin- wave excita- 
tions acquires a gap so that the long-range order is rather 
destroyed by DW proliferation at finite temperatures. To 
fix the ideas, one can think the field fi(x) to be associ- 
ated with localized excitations (DWs) while 4>{x) being 
associated with spin waves. Within a distance separat- 
ing two successive DWs one can assume (n^)^ ~ 1 and 
l^j^np ~ 0. Therefore, up to quadratic terms, the 0-field 
Hamiltonian reads: 



A/2 
-A/2 



dx . 



(15) 



A being the average distance between two successive DWs 
(see Appendix [B] for the derivation of Eq. (IT5|) '). Follow- 
ing the typical procedure of the renormalization group 



(RG), we express the Hamiltonian (fT5)) in the Fourier 
space 



^^^xEi / {'■^l' + D]\MQ)\'dq, (16) 



(17) 



and apply equipartition 



{\Mq)\') = 



Jg2 + 2D 



The integration of the /asf-fiuctuating excitations in the 
range -p- < 9 < f-, with L' = L + dL, yields 



T 



tt/L 



dq 



T 



dL 



(18) 



where ~ J/2D. The prefactors of {dxfi) and (n^) in 
Eqs. (fTH) together with Eq. (fTS|) define the RG equations: 



J(L + dL) = J{L) [1 - {4,1) 
D{L + dL) = D{L) [l-m 



(19) 



As only short-range order is present in a Id system, it is 
reasonable to integrate Eqs. (fT9|) only up to the average 
distance separating two successive DWs, A. The inverse 
of A can be identified with the average density of DWs. 
A low-temperature estimate of the DW density is given 
in Ref. [36|: 



1 1 4£u, ( 



(20) 



where =2%/2Jd again. Note that in Ref. [M the DW 
density is given in terms of bare constants J and D. Our 
goal is to see how A has to be corrected at finite temper- 
ature to account for spin-wave renormalization. This has 
been done with the following procedure: 

1. fix the temperature T 

2. compute numerically the renormalized constants at 
an infinitesimal dL starting from the hare values of 
J and D 

3. compute K(dL) from Eq. H20p using the renormal- 
ized constants J{dL) and D{dL) 

4. iterate the procedure till L + dL > A{L) 

5. at this point we defined A(T) = A(L) and anal- 
ogously the renormalized constants: D{T), J{T), 
e^{T) and w{T). 

Apart from a multiplicative constant, the final value of 
A(T) and the correlation length, f , are supposed to de- 
pend on temperature likewise. 

Before proceeding with the analysis of the RG equations, 
we specify the convention on the notation we use. In this 
section, we introduced a dependence on L in the quanti- 
ties D{L), J{L), £w{L) and w{L) in order to renormalize 
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FIG. 3. Color online. The ratio between renormalized con- 
stant D{L) and the bare constant D{0) is plotted as a func- 
tion of some selected intermediate values of L, for different 
values of D{0) = 0.05 (black squares) ,0.1 (red circles), 0.2 
(blue crosses) and different temperatures ew{0)/T — 0.0625, 
0.106, 0.149, 0.193, 0.236. The renormalization flux stops at 

L = A(r). 



them. This dependence is somewhat technical. The rel- 
evant values of those renormalized quantities are D{T), 
J{T), £w{T) and wlT), namely the values assumed when 
the renormalization stops. These last values only depend 
on temperature, not on L anymore. In the other sec- 
tions, when not specified difi'erently, we always refer to 
the bare constants; these equal the values of the corre- 
sponding renormalized quantities at T = and those of 
the "technical" D{L), J{L), £^{L) and w{L) for L = 0. 



FIG. 4. Color online. Same data as in Fig. [3] but plotted vs. 
L/w{0). 



temperature) and focus on the renormalization of D with 
increasing L. In order to compare the renormalization of 
the anisotropy energy with L for different values of the 
bare D(0), we plot the ratio D{L)/D{0) (Fig. ^. The 
different temperatures have been chosen so that the ra- 
tios ew{0)/T are the same for different initial anisotropy 
values. Fig. m shows how all the curves corresponding to 
the same eu,(0)/r ratio but different D{0) collapse onto 
a single one when plotted as a function of L/w{0) (in- 
stead of L). Such a data collapsing evidences that due 
to the symmetry of Eqs. (|22p the renormalization flux, 
i.e. the relative change of renormalized constants with 
respect to the bare ones, eventually depends only on the 
initial values e^(0)/r and L/w{0). In other words, what 
matters are 



C. Symmetry properties of 
the renormalization flux 

By noting that iwD = one can rewrite the average 
of the 0-field components as: 



T w 



1 



dL 



(21) 



From the above expression, one can immediately see that 
((/)^) depends on L/w{0) and Ew/T only. This symmetry 
property is directly transferred to the RG equations. This 
can be made more transparent by rewriting Eqs. (jl9p as 



dj_ 

J 



T w 1 dL 



dP _ o T w 



1 dL 



(22) 



Note that the renormalized constants will, however, de- 
pend on the values they take at zero temperature J(0) 
and D{0), i.e. the bare constants. Without loss of gen- 
erality, we can set J = 1 (as a unit for the energies and 



• the temperature in units of the DW energy, (0) 

• and length scales in units of the DW width, w{0). 

Besides that, from Eq. (|20p it is clear that A/w only de- 
pends on £w/T . Therefore, applying the same argument 
as for the renormalization flux displayed in Fig. 21 we 
expect A(T)/u;(0) to be a function of e„(0)/T only. The 
correlation length is expected to fulfill the same scaling 
property as A: S^/w{0) has to be a universal function of 
eiii(0)/r. In the following we will see that this property 
holds true provided that £^,(0) is replaced by the DW 
energy obtained from the discrete-lattice calculation £dw 
(we remind that Fig. [1] highlights some deviations of 
from e^u) in the vicinity of £'/J = 2/3). 



D. Transfer Matrix versus 
Polyakov Renormalization 

As stated before, the temperature dependence of the 
correlation length is expected to follow the behavior of 
A(T). In Fig.O we compare ^(T) computed by means of 
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FIG. 5. Color online. Logarithm of the correlation length vs. 
e^{0)/T: TM calculation for D = 0.05 (solid red line); A(r) 
multiplied by a constant to match the TM results at low T 
(circles); analytic result for D = given in Ref. 12 (dashed 
line). The straight lines are guide to the eye to highlight the 
change of slope due to spin-wave renormalization. The dot- 
ted line corresponds to the renormalized DW width ln[u;(r)]. 
Inset: Logarithm of £,/w[0) vs. £d„(0)/r for D = 0.1, 0.2, 
0.3, 0.4, 0.5, 0.6 (same data as in Fig. [5}. 



a universal scaling function of e^(0)/T. However, the 
whole RG approach is based on the continuum approxi- 
mation, Eq. pT| . which is expected to hold for L)<C J, i.e. 
when DWs are significantly broad. But in fact, the dis- 
crepancies between £u)(0) — 2\/2JD and the correspond- 
ing values obtained from the discrete-lattice calculation, 
Edwi are significant for most of the D and J used to pro- 
duce the curves in Fig.[5J For D/J<2/Z but not D<^J, 
we can tentatively extend the validity of the scaling prop- 
erty of ^ beyond the continuum limit by replacing e„ with 
the DW energy obtained from the discrete-lattice calcu- 
lation. Therefore we propose that 




with / universal scaling function. In the inset of Fig. [5] 
we check the validity of Eq. ^ by plotting ^(r)/w;(0) 
as a function of the ratio Sdw/T for J = 1 and different 
values oi D — 0.1 — 0.6. The correlation length data are 
the same as the ones plotted in Fig. [2] The data collaps- 
ing predicted by Eq. ([25)) is indeed observed in the inset 
of Fig. O Remarkably, the log-linear plot of the scaling 
function / {sdw/T) vs. £dw/T does not show a constant 
slope as a consequence of the intrinsic temperature de- 
pendence of characterizing the broad-wall regime. 



the TM technique with A(T) as a function of ew(Q)/T. 
We chose the specific value D = 0.05 because for larger 
values of D the correction to the DW energy due to the 
discreteness of the lattice is not negligible and the identi- 
fication edw — 2\/2DJ is not totally justified (see Fig.[T]). 
The values of A(r) (crosses) have been shifted by a con- 
stant factor to match the TM results (line-symbols) at 
low temperature. Indeed, the temperature behavior of 
A(r) closely follows that of the correlation length till rel- 
atively high temperatures. The change in the slope oc- 
curring at intermediate temperatures is well reproduced 
by the RG calculation. This means that spin- wave renor- 
malization is the main physical reason for the reduction 
of with increasing temperature observed in the broad- 
wall regime (dotted lines in Fig [2]). In the RG language, 
such a change in the slope can also be interpreted as a 
crossover towards the temperature at which D{T) van- 
ishes. This phenomenon is the analogous of the reori- 
entation transition in magnetic films^^. The progressive 
vanishing of the uniaxial anisotropy with increasing tem- 
perature reflects in the fact that for the highest temper- 
atures reported in Fig. [5] the TM calculation recovers the 
analytic solution for the isotropic Heisenberg chain, with 
D{0) — (solid line)^^. In such a region, the behavior 
of A(T) deviates from the one of the correlation length. 
However, at such short distances the RG treatment loses 
its meaning since the length scale at which the renor- 
malization stops, A(T), becomes of the order of the DW 
width 'w{T) (dotted blue line) or even smaller. 
According to the RG analysis discussed in subsec- 
tion |lTrC] one expects that the ratio ^/w{0) be given by 



IV. DYNAMIC PROPERTIES: TQMC 

In order to study the dynamics of the model defined 
by Eq. ([1]), we used a time-quantified Monte-Carlo al- 
gorithm (TQMC)i^ii^ fixing N = 100 and implement- 
ing both periodic and open boundary conditions. In the 
TQMC scheme, Monte-Carlo steps (MCS) are mapped 
into real-time units through the following relationii, 



At[TK] = At[MCS] 



(24) 



This relation gives the variation of real time in units of 
the damping time tk where R is the size of the cone 
used for updating single-spin configuration o^^''*^ . Note 
that the factor that relates MCS with the real time is di- 
vided by the temperature T, therefore the time spanned 
in a given simulation increases when the temperature 
is decreased, provided the factor DR^/{AOT) is small 
enough. We set R = 0.025 in all simulations. Thus, 
for typical values D/J = 0.1 and T/J = 0.1 one has 
DR^/{AQT) ~ 1.5 X 10"^ On the other hand, the damp- 
ing time is given by 



1 



1 



(25) 



where fJ^oHx is the anisotropy field, a the adimensional 
damping constant and 7 the gyromagnetic ratioi^. Since 
the anisotropy field associated with Hamiltonian ([I} can 
be expressed as hoHk = 2D/gfiB and 7 = giis/fi, 
then jfioHK = 2D/h. Henceforth, the quantity h/2D 



TABLE I. Parameters and physical quantities. 



D/J 






r/J range 


S 




simulation 


0.1 


0.894 


2.23 


0.116 


- 0.26 


86 - 


5 


Relaxation 


n 9 
u.z 




i.Oo 


0.153 


- 0.3 


ou - 





Relaxation 


0.3 


1.549 


1.29 


0.181 


- 0.33 


83 - 


5 


Relaxation 


0.1 


0.894 


2.23 


0.01 - 


0.04 


86 - 


5 


Diffusion 


0.2 


1.264 


1.58 


0.014 - 


0.056 


89 - 


5 


Diffusion 


1 


2 


1 


0.1 - 


0.25 


T-IO'^ - 


165 


Diffusion 


2 


4 


1 


0.25 - 


0.48 


536 - 


10 


Diffusion 



The DW width w and the correlation length ^ are given in 
lattice units. The latter has been rounded to integer values. 



will be assumed as time unit. In real SCMs D is of 
the order 10 K. We used an intermediate value for the 
damping constant, a = 0.25, therefore tk is in the 
range of picoseconds {Ti = 7.6 K ps), with these param- 
eters. Moreover, the maximum number of MCS that 
we spanned in our simulations is of the order of 10® 
[MCS], which would correspond to a total time ~ 0.1/is 
for D/J = T/J. For TQMC simulations we set the fol- 
lowing values for the ratio between the anisotropy and 
the exchange strength £'/J = 0.1, 0.2, 0.3, 1 and 2. These 
ratios are comparable to values of real SCMs^i^ and cor- 
respond to different regimes of the DW profil o^^'^^ (see 
Fig. [T|). We studied magnetic relaxation for the selected 
values D/J — 0.1,0.2 and 0.3, all falling in the broad- 
wall regime. The temperatures were ranged between the 
values indicated in Table HI In relaxation simulations, 
such values of T were chosen in order that the corre- 
lation length was smaller than the system size, even if 
for the lowest temperatures ^ and N = 100 are of the 
same order of magnitude. The applied field was cho- 
sen so that H = 0.05Z) in every simulated relaxation. 
In the next section we will present results concerning a 
detailed study of the trajectory displayed by DWs be- 
cause of thermal fluctuations. This investigation was 
also performed by means of TQMC simulations for dif- 
ferent temperatures and anisotropy values. For the dif- 
fusion studies in the broad-wall regime we set D/J~0.1 
and D/J = 0.2; temperatures were chosen in the ranges 
T/J = 0.01 - 0.04 and T/ J = 0.014 - 0.056, respectively 
(see Table |T| . For the studies in the sharp- wall limit we 
used D / J = 1 and 2 while temperatures were chosen in 
the ranges r/J = 0.1-0.25 and T/ J = 0.25 - 0.48 respec- 
tively. Table [1] summarizes the parameters and physical 
quantities related to the different types of simulations. 



A. Relaxation curves 

For studying how the magnetization relaxes towards 
equilibrium we used the following protocol. In each sim- 
ulation, the system was cooled down in zero applied field 
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Time [h/2D] 
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FIG. 6. Relaxation curves at different temperatures for 
D/J = 0.2. Inset: FC curve obtained by TM calculations 
(solid line); open circles are the values of Mo obtained by 
fitting the relaxation curves with Eq. (|26|) at different tem- 
peratures. 



at a cooling rate = 5 • 10^^; this means that the 
temperature was decreased according to the expression 
T = To — rot, with To/ J = 1 being the initial temperature 
and t the time measured in MCS. Once the temperature 
at which we wanted to simulate relaxation was reached, 
we let the system equilibrate in zero field. Further on, 
we switched the field on and let the magnetization evolve 
towards its equilibrium value. The equilibration time in 
zero field was chosen to be equal to the relaxation time 
(in H ^ 0). Fig. [H] shows the relaxation of the magne- 
tization under a homogeneous applied field of strength 
H/J = 0.01 along the z direction for different tempera- 
tures. The anisotropy value is D/J = 0.2. 
All these curves can be fitted well with a stretched expo- 
nential, like in Ref. [2^: 



MS) = Mo 



1 



(26) 



For all the relaxation processes we found that 0.8 < a < 
1. In the inset of Fig.[6]we report field-cooling curve (FC) 
as a function of T/J obtained by TM calculations for an 
infinite system (full continuous line) . The applied field is 
the same as in relaxation simulations, H/J = 0.01. The 
points represent the values of Mq obtained from fitting 
the computed relaxation curves with Eq. (j26l) . The very 
good agreement confirms that the magnetization indeed 
relaxed to its equilibrium value in all the simulated re- 
laxation processes. Fig. [7] shows the relaxation times, r, 
obtained by fitting different relaxation curves, analogous 
to the ones shown in Fig. [6l for three different values of 
the anisotropy D/J = 0.1, 0.2 and 0.3. The relaxation 
times are plotted in a log-linear scale as function of the 
inverse temperature normalized to the corresponding DW 
energy. An Arrhenius dependence of r on the tempera- 
ture is indeed evidenced for any value of D and for both 
open and periodic boundary conditions (o.b.c. and p.b.c. 
respectively). The values of r obtained with o.b.c lie on 
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limit. For the latter case we will provide a numerical con- 



firmation of the phenomenological law Ds 



-D/T 



(see 
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FIG. 7. Color online. Arrhenius plot of the relaxation times 
obtained at different temperatures and for three values of the 
anisotropy, D/J — 0.1 (full circles), 0.2 (squares) and 0.3 
(triangles), with open boundary conditions. The dashed line 
is guide to the eyes and has slope 1.1. Empty circles: periodic 
boundary conditions for D/J = 0.1. 



the same line whose slope is roughly 1.1. The same calcu- 
lation was repeated with p.b.c. for D/J — 0.1 only. The 
corresponding points (empty circles in Fig. [T]) can be as- 
sumed to lie on the same line as for o.b.c. up to the value 
£w/T ~ 6.5. The last point at lower temperature devi- 
ates significantly from the dashed line with slope 1.1. At 
this temperature, the correlation length is comparable to 
the system size, which explains the discrepancy between 
the calculation with p.b.c and o.b.c. (see Table |I]). More 
specifically, it is ^ = = 100 for eJ/T = 7.86. On the 
other hand, for s^/T <6.5 the relaxation time computed 
with open and periodic boundary conditions depends on 
the temperature likewise. This fact suggests that for such 
temperatures the relaxation is not affected by the details 
of boundary conditions, namely it is a bulk process. 
We cannot provide any trivial explanation for the uni- 
versal slope observed for e^/T < 6.5: A,- ~ l-lSw We 
remark that in the sharp-wall limit the formula given in 
Eq. ([5]), At- = 2Aj -I- A^, has been confirmed by a num- 
ber of experiments on SCMs^""— i^^. Due to the fact that 
for sharp DWs Aj = £u, (with £u, = 2J), the same formula 
adapted to the broad- wall limit would predict a slope of 
about 2 + D/ew The value obtained by fitting the data 
in Fig.[7]is therefore about 50% smaller than what would 
predict Eq. ([8]). We will come back to this important 
point when comparing our results with the few available 
experimental data on SCMs in the broad-wall limit. 



V. DOMAIN- WALL DIFFUSION 

In this section we discuss the kind of trajectory dis- 
played by DWs at finite temperature. Particular atten- 
tion will be given to the temperature dependence of the 
diffusion coefficient both in the broad- and sharp-wall 



Eq. ([7])) which - to the best of our knowledge - was still 
missing in SCM literature. The Arrhenius-like depen- 
dence of the diffusion coefficient highlights that each ele- 
mentary move of a DW occurs - in the average - through 
a thermally activated mechanism, in the sharp- wall limit. 
We will show that this is not the case in the broad-wall 
regime. 



A. Analysis of domain-wall trajectories 

As already mentioned, the assumption that DWs per- 
form a random walk induced by thermal fluctuations is 
the basic ingredient to relate the correlation length to re- 
laxation time. Such assumption can be verified directly 
in TQMC by analyzing the microscopic configurations 
explored during a simulation. Details about this anal- 
ysis can be found in Appendix [Cj The outcome is the 
average trajectory followed by each DW at a given tem- 
perature. We found that, the diffusion relation (a;^) cx i is 
not obeyed at short times (see Eq. (|C1[) ). However, DW 
trajectories could be fitted well with a more general ex- 
pression which describes a random walk with correlated 
steps^: 




(27) 



Tc is the characteristic crossover time from the ballistic 
regime at short times to the diffusive regime at longer 
times. When t <ti Tc we are in the regime of correlated 
steps. Expanding Eq. ([27|) accordingly, for ^/tc ^ 1, one 
obtains the ballistic relation between the displacement 
and time: 



When t ^ Tc, the diffusion equation is recovered: 



Tc 



(28) 



(29) 



with Ds being the diffusion coefficient {2Ds = a'^/rc). 
By fitting the mean-square displacement of DWs with 
Eq. ((27| both the diffusion coefficient and the crossover 
time, Tc, can be obtained. In Tablellllwe report the values 
of such parameters for different temperatures, D/J = 0.2 
and H = 0. Note that both Dg and Tc decrease when 
the damping constant increases. For low applied fields, 
like those used in relaxation simulations, Ds and Tc are 
independent of the field itself. 

In conclusion, over some time window - generally 
larger the lower the temperature is - each DW performs 
a ballistic motion before the genuine diffusion process 
starts. 
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TABLE II. DiflFusion parameters. 



T/J 


a 




a 


2Ds 


0.14 


0.25 


85.0 


4.04 


0.44 


0.18 


0.25 


16.7 


2.01 


0.49 


0.22 


0.25 


3.97 


1.10 


0.52 


0.14 


4.00 


32.3 


2.01 


0.35 


0.18 


4.00 


7.00 


1.09 


0.41 


0.22 


4.00 


2.12 


0.66 


0.45 



Temperature dependence of the diffusion 
coefficient 



In order to obtain quantitative results about the 
temperature dependence of diffusion coefficient (see 
Eq. (|29p) we introduced a DW at the center of the spin 
chain with anti-periodic boundary conditions (the spins 
at each end were kept anti-parallel to each other along 
the z axis: — 1 and 5^ = — 1). In each numerical 
experiment, we thermalized the system at a given tem- 
perature and then followed the DW trajectories. The 
regimes with £'/J>2/3 and £'/J<2/3 will be analyzed 
separately. 

Sharp-wall regime - In Fig. [5] we plot in a log-linear 
scale the diffusion coefficient as a function oi D/T for 
anisotropy D/J = l and 2. In this scale the points sim- 
ulated can be fitted well by a straight line indicating an 
Arrhenius behavior, Ds ^ e~'^^l'^ . This fact confirms 
the validity of the expression given in Eq. ([7]) for the dif- 
fusion coefficient in the sharp-wall regime: the slope is 
Aa = 0.891? and Aa = 0.971) for D/J=l and 2, respec- 
tively. The slight reduction of the energy barrier with 
respect to the prediction of Eq. ([7]) (A^ = £)) may be due 
to the high-temperature points in Fig. [5] Even for the 
attempt frequency of a single nanoparticle with uniaxial 
anisotropy^, an Arrhenius behavior is expected only in 
the limit D/T <^ 1. Another possibility is that A^ be- 
comes smaller than D as the crossover ratio, Z3/J = 2/3, 
is approached. In this sense one would expect the renor- 
malizing effect of spin waves to be more important for 
D/J = 1 than for D/J = 2 (see Sec. Hill. Therefore, 
within the numerical accuracy, our TQMC simulations 
confirm the phenomenological law proposed in Ref. .25i 
and [39! with IS.a^D. Broad-wall regime - The tempera- 
ture dependence of the diffusion coefficient changes when 
the anisotropy-to-exchange ratio is reduced. In Fig.lHlwe 
plot the diffusion coefficient vs. temperature for the 
two ratios £'/J = 0.1 and 0.2 (broad-wall regime). The 
temperature is now expressed in units of the DW energy 
given in Eq. In this case the diffusion coeffi- 

cients increase linearly with temperature. This behavior 
is at odds with the Arrhenius dependence predicted by 
Eq. ([7]) for sharp DWs: it rather reminds the behavior 
of the diffusion coefficient of a massive particle in a vis- 
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FIG. 8. Color online. Temperature dependence of the diffu- 
sion coefficient for D=l and 2 (sharp- wall regime). 
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FIG. 9. Color online. Temperature dependence of the diffu- 
sion coefficient for two values of the anisotropy D = 0.1 and 
0.2 (broad- wall regime). The inset shows the scaled curves. 



cons medium. The slope of the diffusion coefhcient as 
a function of T /e^ is larger the smaller the anisotropy 
is. Following analogous considerations to those that al- 
lowed us to derive the scaling relation in Eq. (|23p . we 
can attempt to propose a scaling ansatz for the diffusion 
coefficient. Note that the units of Ds are square unit 
lengths divided by a unit time. The time unit assumed 
throughout the paper is (7oMo^fe)~^ = h/2D. This unit 
corresponds, in general, to different physical time scales 
for different values of D. Nevertheless, this specific de- 
pendence on D has been already eliminated from Ds by 
measuring the time in units h/2D. Then we need to ex- 
press the lengths in unit of w and the energies in units 
Sdw, which yields 



DsiT) 



(30) 



From Fig. [9] and from the analogy with a particle in vis- 
cous medium, we conclude that the scaling function g is 
just a line passing through the origin in the T-Ds plane 
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so that 

5# = A^. (31) 

A being some constant with units 2D/fi. In the inset 
of Fig. [5] we plot Ds/uP' vs. T/sdw using the same data 
as in the main frame. Indeed, the scaling prediction of 
Eq. ^ is well-obeyed, giving A = 0.17 [2D/h]. 

VI. PHENOMENOLOGICAL ARGUMENTS 

The analysis performed for both static and dynamic 
properties allows stating that the relevant energy scales 
in our problem are the DW energy and the width 
and position of the spin-wave spectrum. Such energies, 
and their relationship with the thermal energy, determine 
the physics of SCMs described by the model Hamilto- 
nian ([T]). The full spin- wave spectrum can be obtained 
by linearizing the Landau-Lifshitz equation correspond- 
ing to Hamiltonian ([1]) . The energy of a spin wave with 
frequency cj and wave vector q is 

huj{q) ^ [2 J (1 - cos(g)) + 2D] . (32) 

At low enough temperature, in a spatial region delimited 
by two DWs one essentially has | S*^ | = 1 . In the left panel 
of Fig. [10] the spectrum of fluctuations, huj{q), is plotted 
vs. q ioi D — 2.5 J (the sharp- wall limit). The dashed 
horizontal line indicates the corresponding DW energy 
£w = 2 J. Clearly the energies of the two family of exci- 
tations - spin waves and DWs - are well separated from 
each other. Moreover, the genuine Id character of a spin 
chain is evident when the energy of thermal fluctuations 
is lower than the DW energy: T < e^,. At these tem- 
peratures some short-range correlations develop, i.e., ^ 
exceeds some lattice units. The dotted horizontal line 
highlights a realistic reference for such a thermal en- 
ergy. On the right panel, the same plot is displayed for 
D = 0.2 J corresponding to broad DWs. In this case the 
dashed horizontal line, representing the DW energy e^,. 




q 



FIG. 10. Color online. The solid lines represent the energy 
of spin waves fiu}[q) (in J units) as a function of the wave 
vector q. Horizontal lines correspond to the energy of one 
DW in the sharp- and in the broad-wall regime: for the left 
panel D/ J = 2 while for right panel D/ J — 0.2. The horizontal 
dashed line represents an indicative temperature smaller than 
Em in both cases (see the text). 



passes through the spectrum of spin waves. As a con- 
sequence, in the broad-wall regime one expects the in- 
terplay between DWs and spin waves to affect crucially 
the finite-temperature properties of the system. On the 
other hand, as in Fig. (TU] (left) the spin-wave spectrum 
lies well above and T, spin waves are expected to play 
no essential role in sharp- wall limit (when T <Sw <2D). 
The fact that computed in Sect. Hill is independent of 
T for sharp DWs while it effectively depends on tempera- 
ture for broad DWs, confirms the heuristic argument evi- 
denced by Fig. 1101 In particular, we have shown through 
Polyakov renormalization that the interplay between spin 
waves and DWs at finite temperature gives a quantita- 
tive explanation for the dependence of on T in the 
broad- wall limit. Parenthetically, we note that in the 
time domain it is easy to see that averaging over spin 
waves (the scalar fields 0q in the language of Polyakov 
renormalization) corresponds to an integration over fast 
fluctuations. In fact, according to Eq. ([5^ the typical 
time periods of spin waves, 27r/w(g), are of the order of 
our reference time unit U/2D or even smaller. The typ- 
ical time scale for the creation or annihilation of DWs 
is, instead, of the order of the relaxation time r, i.e., 
several orders of magnitude larger than h/2D. Thus, in 
the experimental situation relevant for SCMs it is clearly 
T'^l/uj(q), meaning that n{x) represents a slow varying 
field and (t>{x) a fast varying field. 

For what concerns the relaxation time we cannot pro- 
vide an effective argument as Polyakov renormalization 
to justify our numerical findings. However, it seems rea- 
sonable that the very same interplay between DWs and 
spin waves affects the temperature dependence of the re- 
laxation time in a similar way as it affects the correlation 
length. A natural consequence of this is that Eq. ^ does 
not necessarily hold true in the broad-wall limit. One has 
to be very cautious even in trying to generalize the rela- 
tion A^ = 2A5 -I- A^ (Eq. (O) to the case of broad DWs. 
In this regard, we recall that for D/J<2/3 

1. A J depends on the temperature range in which it 
is measured 

2. the time window in which the DW motion is ballis- 
tic, and not diffusive, becomes larger while lowering 
the temperature 

3. in the temperature range that we investigated with 
TQMC it is Aa = 0, meaning that each single DW 
move is not thermally activated. 

For < T < 9 Ew and £^ < N, the numerical data re- 
ported in Fig. [7| suggest that Eq. ([5]) has to be modified 
into 

A, = l.le^ (33) 

for broad DWs. In order to check this prediction against 
experiment we refer to two Mn(III)-based spin chains 
with D/J<2/3 (broad DW) reported in Ref . 'sj and IH. 
Both chains are better described by the Seiden model^ 
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with anisotropy (on the classical spin) rather than the 
Heisenberg model in Eq. ([T]). In fact, the Mn centers 
(5 — 2) alternate with an organic radical (TCNE or 
TCNQ) whose magnetic contribution is essentially the 
same as a free electron: spin s = 1/2 and Lande factor 
gs =2. In the broad- wall limit, the Seiden model with 
anisotropy can be mapped into the Heisenberg model de- 
scribed by Hamiltonian ([1]) with an halved exchange cou- 
pling^. For the Mn(III)-TCNE spin chain^l, we have 
estimated in Sect. IIIII g„, ~ 100 K. Adapting to our con- 
vention the values of D and J given in Ref. IH we ob- 
tain ~ 86 K for the Mn(III)-TCNQ spin chain. Thus 
Eq. ([33)l would predict the following activation barriers 
for relaxation 

• for the Mn(III)-TCNE spin chain ~ 110 K, to 
compare with the experimental value A^^*' = 117 K 

• for the Mn(III)-TCNQ spin chain A^ ~ 95 K, to 
compare with the experimental value A^^p = 94 K. 

The prediction of Eq. p3p agrees well with the measured 
energies in both cases. We note, however, that in exper- 
iments the validity of the empirical formula A,- = l.le„ 
seems to extend down to a temperature region in which 
£^ > N. In fact, the lowest temperature for which the 
A«^P ^ 94 K for the Mn(III)-TCNQ spin chaiii^a is 
roughly T = 4.5 K. The corresponding correlation length, 
extrapolated from Fig. ^ should be of the order of 10^ 
Mn(III)-TCNQ units. The same estimate gives a corre- 
lation length of the order of 10^^ units for the Mn(III)- 
TCNE spin chain^^ As already stated in Sect. [TTTl de- 
fects and dislocations typically limit the length of spin 
chains to 10^ — 10^ units (see the dotted horizontal lines 
in Fig.[2])- Therefore, at the lowest temperatures at which 
Eq. ([55)1 seems to apply, ^ should exceed the average dis- 
tance between two defects in both molecular spin chains. 
A conclusive analysis would require a more accurate fit- 
ting of the model parameters to the experimental data 
for each sample. At this stage, we have no qualitative 
explanation nor a numerical confirmation for the validity 
of Eq. ([55]) in the regime ^ > A^. Simulating a relaxation 
experiment at lower temperatures, in the region where 
^ ^ N, is computationally very expensive due to the 
Arrhenius dependence of the relaxation time. This issue, 
indeed, deserves further theoretical investigation but this 
is beyond the scope of the present work. 

VII. CONCLUSIONS 

We studied a model paradigmatic for classical spin 
chain or magnetic nanowires with uniaxial anisotropy 
and identified two distinct regimes for static and dynamic 
properties. Such difl:erences in the finite-temperature be- 
havior are closely related to the thickness of DWs at zero 
temperature. We distinguished, accordingly, between the 
sharp- and broad-DW regimes. In the sharp-wall regime 
{D/J>2/3) the correlation length obtained by TM cal- 
culations shows an activated behavior as function of the 



inverse of the temperature. The corresponding activa- 
tion energy is equal to the DW energy A^ —e^ {£w = 2 J 
in this regime). In fact, for large anisotropy-to-exchange 
ratios, DWs extend only over one lattice spacing so that 
the anisotropy energy does not affect two-spin correla- 
tions. At variance, when DWs develop over several lat- 
tice units (Z?/J<2/3), the correlation length still shows 
an exponential divergence with the inverse of the tem- 
perature, but with a temperature-dependent A^ . At low 
temperatures the activation energy is larger A^ = O.Qew 
(here e^u = 2\/2_D J), whereas at higher temperatures 
it is A^ = 0.6e^. The first result agrees with a low- 
temperature expansion available in the literature^. Be- 
sides that, we provided a physical argument - based on 
Polyakov renormalization - to justify the 30% of reduc- 
tion of A^ observed at high temperature. This allowed 
us to conclude that the appearance of the lower activa- 
tion energy at higher temperatures is due to the interplay 
between DWs and spin waves. This interplay is not signif- 
icant in the sharp- wall regime where static properties are 
practically determined by the energy cost to create a DW 
at zero temperature. The reason why physics is remark- 
ably different in the sharp- and broad- wall regime lies on 
the relative difference among the energy scales involved 
in the problem and the thermal energy corresponding to 
temperatures at which short-range correlations develop 
(see Sect. EH. 

In SCMs, relaxation is usually assumed to be driven by 
DWs diffusionSi"— Based on this assumption, the acti- 
vation energy for the correlation length, Aj, and that of 
the relaxation time, A^, are then related with each other. 
We tested, with TQMC simulations, that DWs indeed 
perform a random walk for time intervals much longer 
than the typical precessional time of a single spin and de- 
termined the temperature dependence of the diffusion co- 
efficient. In the sharp- wall regime, we found that the dif- 
fusion coefficient Dg follows an Arrhenius behavior with 
an activation energy close to the anisotropy value, D, as 
assumed in most of the experimental works^"— i^^i^. In 
the broad- wall limit, the diffusion coefficient does not fol- 
low an activated mechanism but it rather grows linearly 
with the temperature, reminding the behavior of a par- 
ticle in a viscous medium. The results of this analysis 
confirm the robustness of the random- walk argument re- 
lating the correlation length to the relaxation tim o^^'^^" — 
and suggest a dynamic critical exponent z = 2. Neverthe- 
less, in the broad- wall regime, the relation between the 
A^ and A^ is not trivial due to spin-wave renormaliza- 
tion. As a consequence, the joint theoretical and exper- 
imental characterization of SCMs falling in this regime 
cannot be based on the simple Glauber model^S or gen- 
eralizations of itii24-^i39^ 

The symmetry of renormalization-group equations sug- 
gests the existence of scaling laws specific to the broad- 
wall regime: the na tural unit for length scales is the DW 
width w = -y/ J/ 2D while for energies it is = 2V2DJ 
(the DW energy at zero temperature). Thus, one ex- 
pects that the physical observables which depend only 
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on these quantities be universal functions of properly 
rescaled variables (e.g. e^jT). Static TM calculations 
and dynamic TQMC simulations confirm the validity of 
this conjecture for the correlation length, the diffusion co- 
efficient of DW motion and the relaxation time (for the 
latter, scaling is obeyed apart from a residual tempera- 
ture dependence in the time unit intrinsic of the TQMC 
method). 

In view of possible magneto-storage applications, in- 
creasing the thermal stability of the SCMs would be de- 
sirable. To this aim. we note that designing novel com- 
pounds with a larger J as possible would not be a good 
strategy for two reasons: for Z?/ J< 2/3 i) the DW energy 
is the sole quantity which controls thermal stability and 
it scales as ^ ^fj, instead of ^ J like for sharp DWs; 
ii) spin- waves renormalization progressively lowers the ef- 
fective energy barrier for relaxation as the temperature 
is increased (by renormalizing e„). 

Temperature is often neglected in models employed 
to study the current-induced DW motion in magnetic 
nanowires^"— or it is taken into account in the phe- 
nomenological parameters of the Landau-Lifshitz-Gilbert 
equation^. The basic hypothesis is that the considered 
nanowire behaves as a 3d magnet below its critical tem- 
perature^i^. However, recent experiments^""— suggest 
that Joule heating may induce the formation of domains 
in the nanowires, which highlights the restorations of a 
genuine Id magnetic character. In this situation, the DW 
trajectory may result from a delicate combination of ther- 
mal diffusion (stochastic) and the deterministic motion 
induced by the electric current. Our study of the tem- 
perature dependence of the diffusion coefficient can be 
considered a preliminary contribution to this problem, 
of technological relevanc o^'^^ , that indeed deserves fur- 
ther investigation. Typically in metallic nanowires (Co, 
Ni, Fe, Permalloy) D ~ 1 - 10 K (- 0.1 - 1 meV) and 
J ~ 100 - 500 K (- 10 - 50 meV), meaning that they 
generally fall in the broad-wall regime where the tem- 
perature dependence of any observable is expected to be 
affected by the non-trivial interplay between DWs and 
spin waves. 
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Appendix A: The transfer-matrix approach 

Given a general classical spin-chain Hamiltonian with 
nearest-neighbor interactions 

N 

n = -Y,v{s,, s,+i) (Al) 
1=1 

the partition function Z is given by 

Z = / df^i / df]2 . . . / ef'^^^^'^-^ ePV(s,s,) 

^ _ g/3y(s„,5i)^f^^ (A2) 

where the integrals, dfJi, are performed over any possible 
direction of the unit vectors St and j3 — 1/T. Defining 
the transfer kernel /C®° as 

JC{S,J,+i)^e^^^^-^^+'^ (A3) 

and taking periodic boundary conditions (A^+1 = 1), the 
partition function Z takes the form of the trace of A^-th 
power of }C{Si, Si+i): 

Z = /dr!i/rfr!2...//C(5i, 52)/C(^2, ^3) 

... JC{Sn, Si)dVlN ^Ty{K^] . (A4) 

The computation of such a trace, as well as other physical 
observables, is simplified if one first solves the following 
integral eigenvalue problem: 

J IC{Si, Si+i)'l}jniSi+l)dni+i = \mi'm.{Si) (A5) 

The eigenfunctions fulfill the properties 

ip-m{Si)tpm{Sj) = S{Si - Sj) completeness(A6) 

m 

and 

J ij}n{S)i^m{S) dVl = 5n,m orthonormality( A7) 

where 5{Si — Sj) is the Dirac (5-function and 5n,m 
is the Kronecker symbol. For non-symmetric ker- 
nels, lC{Si, Si+i) 7^ K{Si+i, Si), properties similar to 
Eqs. and hold for the left and right eigen- 

functions, but this is not our case. Using Eq. (|A6[) the 
kernel can be rewritten as 

IC{Si, Si+i) = ^ A„i Tp„iiSt)'ip„iiSi+i) . (A8) 
Combining Eqs. and (K7}i we get 

m 

The eigenvalues A„i are all real and positive, as the ker- 
nel operator (|A3p is a positive defined function of Si and 
Si+i. For symmetric kernels, /C(S'i, Si+i) — IC{Si+i, Si), 
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the reality of the eigenvalues descends from the analo- 
gous of the spectral theorem for real symmetric matri- 
ces. Moreover, it is possible to show that the spectrum 
of (|A5|) is upper bounded so that the eigenvalues A„i can 
be ordered from the largest to the smallest one: 

Ao > Ai > A2 > . . . 

In the thermodynamic limit the asymptotic behavior of 
the partition function (jA9|) is dominated by the largest 
eigenvalue Aq, yielding 



N- 



X 



N 



(AlO) 



Once the largest eigenvalue Aq is known, the free energy 
and its derivative can be computed from the relation F = 
-T InZ. 



1. Pair-spin correlations 

Here we recalli^^ how pair-spin correlations can be 
evaluated by means of the eigenvalues and the eigenfunc- 
tions defined in Eq. (IA5[) . Consider the ^ component of 
the spin located at site i and the v component of the spin 
located at the site (i-l-r), with /x, v — x, y, z. The average 
we have to evaluate is 

...IC{S,-i, S^)S^ICiS^, S^+i)... (All) 

Following the procedure of the previous section we obtain 

J V'mi+,_i (Si+r) Sf+r -0771, + , {Si+r)dfli+r 
^rrii^r-^^i + r+l ' ■ • ^niN -I ,mN ^niN .nil ■ 

Considering all the repeated indices in the Kronecker 
symbols, we have 

Kai I i'm, (St+r) S'f+r "0"!.+, {Si+r)dfli+r ■ (A13) 

If we substitute to Z its asymptotic expansion, 
Eq. (|A10|) . we need to evaluate the products 



^ An / V A„, , , „ / 



(A14) 



for TV — >■ CXI . It is straightforward to conclude that only 
the terms for which rrii^r = will not vanish. Thus, 
pair-spin correlations are given by 

(5f5.V) = E™. !Ms^)s>:i:,nAs^)dn, 

J 1prndSi+r) S';'_f_^1po{Si+r)dflt+r ■ (Al5) 



In the previous formula, i is a dummy index but the 
result obviously depends on the separation between the 
two considered spins, r. Note that not only the Aq ipo but 
all the eigenvalues and eigenfunctions enter Eq. (jAlSp . 
Eq. (|A15|) can be further simplified when the corre- 
lation function is considered 

(sr^r+,,) - E™. {^y\jMSi)st4'mM)dn,f . 

(A16) 



2. Discretization of the unitary sphere 



The eigenvalue problem defined in Eq. (jA5l) can be 
mapped into a linear algebra problem by sampling the 
unitary sphere with a finite number of points. Given a 
generic function of two angles 9 and 0, say f{9, (j>), the 
integral over the solid angle {dft — d(p sin Odd) can be 
approximated as: 



/(</), 9) dn 



p 

E 

h=l 



Whfiuh) , 



(A17) 



where ut represent the special points that sample the uni- 
tary sphere, Wh are the relative weights and P is number 
of points themselves&l. After discretizing the kernel JC 
in this way, the eigenvalue problem \Ab\ transforms into 
the following linear algebra problem 



p 

E 



WhIC{ui, Uh) 1pm{uh) = Xmlpmiui) (Al8) 



which is usually symmetrized as 



K, 



Lh 



y/wiWhlC{ui, Uh) 

^/wK->Pniuh) 



to yield 



p 

h=l 



(A19) 



The number of special points, P, defines the size of the 
matrix that has to be diagonalized. For the calculations 
reported in this work we have chosen P = 72, 128, 200 
(Mc Laren formula 14-th degree^^ and Gauss spherical 
product formulae of 15-th and 16-th degree^^). The com- 
parison between the results obtained for different sam- 
pling allows estimating the accuracy of the calculation 
at a given temperature. We used the routine DSPEV 
of the LAPACK package to solve the eigenvalue problem 
defined in Eq. (jJTO)) . 



Appendix B: Polyakov renormalization 

We rewrite here for convenience the continuum Hamil- 
tonian given in Eq. pip : 



n ^ -JN 



^-\dJ\''-D{S\x)f 



dx (Bl) 
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and represent S{x) as a superposition of tow fields: 

S{x) = n{x) v/1 - 02 (x) + ${x) . (B2) 

As mentioned in the main text, requiring that = 1 

one necessarily has n{x) ■ 4>{x) = so that 4){x) can be 
expressed on a basis orthonormal to n{x): 



c) = ^ ^aCa with a = 1, 2 



(B3) 



with |ea(a;)| = 1—. Moreover, the fact that \n{x)\ = 1 
implies that dxn{x) is orthogonal to n(x) itself, thus 

dxn{x) ^^CaCa . (B4) 

a 

The gradient term in Hamiltonian (|Blll reads 



(B5) 



Exploiting the orthogonality of the basis (n, Ca), the 
derivatives dxfia can be written as 



dx^a = -Can + fabeb 



(B6) 



where fab is an antisymmetric two-by-two tensor whose 
components can be made negligible with a proper choice 
of the reference frame (ei, 62). Finally, one gets 

\dxS\'' = (1 - (^dxfif + i^xCj^-f + E ^aCbCl)aCl>b 

a ab 

a a 

(dxVi^y . 

(B7) 



+ 



The two terms at the second line of Eq. (jB7P vanish after 
spatial integration; while the last term (third line) only 
contains powers of (jja higher than the square after the 
derivation (thus we neglect it). {dx'Paf' represents 
the "kinetic" term of the field (t>{x). Our goal is to re- 
tain all the other quadratic terms in the field 4'{x) and 
use them to perform thermal averages. After thermal 
averaging, (...), we will have 

{Mb) - 5ab{<Pl) (B8) 
which allows rewriting 

(E CaC60a06) = {4>l) J] ^ (0^) {^Oxuf . (B9) 



Besides this, one has 

(0') = (0?+0^>=2(02) 



(BIO) 



Thus, the Polyakov's result is readily recovered 

{\dxS\^) = [l- (02)] (^dxuf + ^((a.0a)') . (Bll) 

a 

Let us consider now the anisotropy term in Hamilto- 
nian (IbTI): 



= {n^-f (1 -cf')+Y, 4>a<l>belet + O (0,,) , 

ab 

(B12) 

where O (0a) stands for terms linear in 0^ that vanish 
after spatial integration; while after thermal averaging 
one has 



•a) 



In this case, we obtain 

((5^)^) = [1-3(0^)] (n^)V(0^) 



(B13) 



(B14) 



Now we want to derive the quadratic Hamiltonian for 
the field written in Eq. (jl5p and used to perform ther- 
mal averages. We first rewrite the Hamiltonian (jBip in 
terms of the fields (j){x) and nix). In doing this, we will 
keep the approximations we made previously: we will ne- 
glect the terms of order higher than the second in 0a as 
well as what is supposed to vanish after the spatial inte- 
gration. A delicate point of the latter approximation is 
the assumption 



Oa<Pb 



= 5ab 



(B15) 



based on symmetry reasons. The simplified continuum 
Hamiltonian reads 



H 



-JN 



J 



\dxn'\^dx 



E 



D 



J 



{dx(t>a) dx 



M'(i-3*i±*i) + ^i. 



dx 



= -JN+^ J \dxn\^dx-D J {n')^dx 



(B16) 



E"^ / {dx4>afdx 



E 



'dx . 



Within a distance separating two successive DWs we as- 
sume 



{n'^Y ~ 1 
I9,n|2~0. 



(B17) 
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which allows expressing the 0-field Hamiltonian as 



100 



Hi 



N 
1 



E 



A/2 



-A/2 



dx, (B18) 



A being the average distance between two successive 
DWs. 



Appendix C: Domain-wall motion and random walk 

In order to check how the choice of boundary condi- 
tions could affect DW diffusion, we performed a prelimi- 
nary analysis of a standard one-dimensional random walk 
constrained to a finite segment. Such a geometrical con- 
straint is equivalent to the one experienced by a DW 
which diffuses in the Heisenberg spin chain with o.b.c. 
In a random walk process the relation between the mean 
square displacement (x^) and the time t elapsed during 
the walk is: 



Tr 



(CI) 



where a is the mean square displacement and Tj, is the 
mean-time for each step. We considered a discrete ran- 
dom walk in a linear system of TV sites using different 
boundary and initial conditions. Since we set the lattice 
unit a — 1 and take Tc = 1 then = 1. In Fig. 111! we 
plot the square root of the mean-square displacement in 
a system of the same size of our spin chain, = 100. We 
used both reflecting and absorbing boundary conditions. 
In addition, we chose for each case two initial conditions: 
in one we started the random walk from the middle of 
the chain and in the other one at a random position. In 
the four studied cases we observed that any significant 
deviation of the relation given by Eq. (|C1|) started at 
(x^) > 10. Below this value boundary conditions are not 
expected to affect the analysis of DW trajectories either. 

For what concerns the spin chain, we began by study- 
ing the diffusion of DWs which were naturally present in 
the system at the chosen temperature (if ^ < N). The 
setup was the same as the one used to obtain the curves 
of Fig. [6l We chose three representative temperatures, 
r/J = 0.22, r/J = 0.18, and T/J = 0.14 using two values 
for the damping constant: a — 0.25 and a — A. Fig. [T^ 
shows the mean-square DW displacement when the sys- 
tem was equilibrated at T/J = 0.14, for D/J = 0.2. The 
data reported in the two figures on the left-hand-side cor- 
respond to zero applied field, whereas for the two figures 
on the right-hand-side the field was iJ/J = 0.01. Instead, 
figures in the first (Figs. [12] (a) and (b)) and in the second 
raw (Figs.[T3(c) and (d)) correspond to the two different 
values of the damping constant. For all the curves, the 
diffusion relation (x^) (x t is not obeyed at short times 
(see Eq. (jCip ). Thus we fitted those trajectories with a 
more general expression which describes a random walk 




100 
Time [Steps] 

FIG. 11. Color online. Square root of the mean square devi- 
ation of a random walker obtained using different boundary 
and initial conditions. 




Time[ll/2D] 

FIG. 12. Average dispersion of the DW trajectories at T/J = 
0.14 for two damping constant and applied fields In Figs, (a) 
and (c) the applied field is H — and the damping constant 
a — 0.25 and a = 4, respectively. In Figs, (b) and (d) the 
applied field H = 0.01 and the damping constant a = 0.25 and 
a = A, respectively. 



with correlated steps^: 



t e'~c - 1 



(C2) 



the meaning of the parameters Eq. (|C2[) is explained in 
the main text in Sect. fVl 

Fig. [13] displays some individual trajectories that were 
used to produce the curves in Fig. [12] Some correla- 
tion emerges for small displacements, of the order of the 
DW width, which confirms the occurrence of a ballistic 
regime for short times. At very short time scales there 
is an uncertainty in the DW position of the order of a 
lattice parameter. In fact, fast fiuctuations of the DW 
structure produce an offset in the initial position of the 
diffusing DW, intrinsic to the method used to detect such 
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a position. 



1500 2000 

Time[ti/2D] 



FIG. 13. Color online. Domain-wall trajectories at r/J = 0. 14 
for two damping constant and applied fields. In Figs, (a) 
and (c) the applied field is H — and the damping constant 
a = 0.25 and a = 4, respectively. In Figs, (b) and (d) the 
applied field ^ = 0.01 and the damping constant a = 0.25 and 
a = 4, respectively. 
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